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We employ the principle of dynamic equivalence between soft-sphere and hard-sphere fluids [Phys. 
Rev. E 68, 011405 (2003)] to describe the interplay of the effects of varying the density n, the tem- 
perature T, and the softness (characterized by a softness parameter v~^) on the dynamics of glass- 
forming soft-sphere liquids in terms of simple scaling rules. The main prediction is the existence of 
a dynamic universality class associated with the hard sphere fluid, constituted by the soft-sphere 
systems whose dynamic parameters, such as the Q-relaxation time and the long-time self-diffusion 
coefflcient, depend on n, T, and v only through the reduced density n* = na%s{n,T, u), where the 
effective hard-sphere diameter aHs(n,T,v) is determined by the Andersen- Weeks-Chandler condi- 
tion for soft-sphere-hard-sphere structural equivalence. A number of scaling properties observed 
in recent experiments and simulations involving glass-forming fluids with repulsive short range in- 
teractions are found to be a direct manifestation of this general dynamic equivalence principle. 
The self-consistent generalized Langevin equation (SCGLE) theory of colloid dynamics is shown to 
accurately capture these scaling rules 
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The formation of colloidal glasses and gels has been 
the subject of intense study during the last two decades 
[l[ . On the other hand, within their overwhelmingly rich 
phenomenology, molecular glass-forming liquids exhibits 
intriguing universal features 0. It has been natural to 
expect that the phenomenology of both, the glass transi- 
tion in "thermally-driven" molecular glass formers, and 
the dynamic arrest transition in "density-driven" hard- 
sphere colloidal systems, share a common underlying uni- 
versal origin [j?]. In fact, interesting scalings of the equi- 
librium dynamics of simple models of soft-sphere glass 
formers, exposed by systematic computer simulations 
0i 0| I provide an initial clue to one possible physical ori- 
gin of this universality, whose fundamental understand- 
ing, however, still constitutes an important theoretical 
challenge. This challenge provides the main motivation 
of the present work. 

An essential aspect of the equilibrium perturbation 
theory of liquids is the static structural equivalence prin- 
ciple This principle states that a fluid at number 
concentration n and temperature T, whose particles in- 
teract through a moderately soft and purely repulsive 
potential u''-^\r) (where is some measure of soft- 
ness), is structurally equivalent to a hard-sphere (HS) 
system with an effective HS diameter uhs — <^Hs{n, T, ly) 
and an effective volume fraction (pns — 4>hs{^iT,v) — 
■Kna^g{n,T,v) This means that the static struc- 
ture factor (SSF) S{k;n,T,v) of the soft-sphere sys- 
tem is given by S{k;n,T,v) w SHsiko'HSi 4>Hs), where 
Snsikcr; <p) is the SSF of the fluid of hard spheres of di- 



ameter a and volume fraction (p. This static structural 
equivalence automatically implies the universality of the 
thermodynamic properties of the class of soft-sphere flu- 
ids defined, precisely, by this iso-structurality condition. 

The dynamic extension of this soft -hard equivalence 
was proposed more recently [9], thus extending the re- 
ferred thermodynamic universality to the dynamic do- 
main, described by properties such as the self interme- 
diate scattering function (self-ISF) Fs{k,t) or the long- 
time self-diffusion coefficient D^. Some implications of 
this universality, on the dynamic arrest scenario of soft- 
sphere systems, have also been discussed [l3|- In these 
discussions, however, temperature was considered con- 
stant, and hence, its role was never emphasized. 

The main purpose of the present work is to demon- 
strate that the same dynamic equivalence principle be- 
comes a much deeper and more powerful fundamental 
tool when conceived as a general dynamic scaling prin- 
ciple in the density-temperature-softness state space of 
these glass-forming soft-sphere liquids. This defines what 
we refer to as the hard-sphere dynamic universality class 
(HS-DUC) and explains, in particular, some of the in- 
triguing scalings observed in the recent simulations on 
model glass- forming soft-sphere liquids 0, 0] ■ 

Let us first refresh the concept of static structural 
equivalence, now in terms of the radial distribution func- 
tion (RDF) g{r; n, T, ly) of a given soft-sphere model sys- 
tem. The fundamental physical notion is that this sys- 
tem behaves essentially as a hard-sphere system in the 
sense that g{r;n,T,v) « gHs{r / (Jhs]4'hs) where 
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9Hs{r/(^',(l>) is the RDF of the HS system. This iso- 
structurahty condition allows one to write the equilib- 
rium thermodynamic properties of the soft sphere sys- 
tem, such as the equation of state p = p(n, T, v), in terms 
of the corresponding properties of the hard sphere fluid. 
For example, the pressure p{n, T, v) can be written, us- 
ing the virial equation of state @], as p(n, T, v)/nkBT = 
Z(n,T,v) « l+^HsgHs{'^+]4>Hs) = Zhs{4>hs)- Using 
the Verlet-Weis prescription to approximate the contact 
value gHs{^-\-] 4>hs) turns out to be equivalent to approx- 
imating the hard sphere compressibility factor Zhs{4>hs) 
by the Carnahan-Starling equation 0, thus finally lead- 
ing to the following approximate but universal mechani- 
cal equation of state of the soft sphere system, 
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(1 - 4>Hsy' 



(1) 



with (pHS — 4>Hs{n,T,v) determined by the iso- 
structurality condition. This universality is nicely illus- 
trated in the simulation results of Ref. Q (see inset of 
Fig. 4 below). 

The dynamic extension of this soft-hard equivalence 
was discussed in Refs. in the context of the dynam- 

ics of colloidal liquids, in which a short-time self-diffusion 
coefficient describes the diffusion of the colloidal par- 
ticles "between collisions" . It is summarized by the state- 
ment that the self-ISF Fs{k,t;n,T, i/) of the fluid with 
soft repulsive potential u^'^^r) can be approximated by 



Fs{k,t;n,T, v) 



Ff'\kaHS.Dh/als-AHs), (2) 



where Ff^\ka, D^t/cr^;<j)) is the self-ISF of the fluid of 
hard spheres of diameter cr, volume fraction cj), and (for 
simplicity) same short-time self-diffusion coefficient D'^ 
as the soft-sphere fluid. As a direct consequence of this 
dynamic universality, it follows that the long-time self- 
diffusion coefficient ^^(ri., T, z^) (= limf^oo((Ar(t))2)/6t) 
of the soft-sphere liquid, normalized as D* (n,T^v) = 
DL{n,T,v)/D°, is given by 



D* (n, T, v) w D*Hs [(j)HS {n, T, v)] , 



(3) 



where D^glc/)] is the corresponding property of the HS 
system. Similarly, let us define the a-relaxation time 
Ta{k; n, T, v) by the condition Fs{k, Tq.) — 1/e, which we 
normalize as t* {ka; 4>,T* , i/) = k'^D'^Ta{k;n,T,iy). The 
dynamic equivalence principle above then implies that 



with 



HS 



T*{k;n,T, v) « T*fjg[kaHS]4>Hs] 
[ka] (j)] refering to the HS system. 



(4) 



Some consequences of the universality summarized 
by Eq. ([2]) were illustrated in Refs. [§] and 
[lOj in the context of the truncated Lennard- Jones 
(TLJ) pair potential with tunable softness, u^'^\r) = 
e [{(Tlrf - 2{a/rY + l] e{(j~r) (with 0{x) being Heav- 
iside's step function), whose state space is spanned by the 



volume fraction (f) = nna^ /6 and dimensionless temper- 
ature T* = ksT /e. These references, however, discussed 
in detail only the limit of moderate softness {v ^ 1), 
in which the strong similarity with the HS potential 
leads to the additional simplification that ansin^T^v) 
becomes n- independent, and given by the "blip func- 
tion" approximation [lo| . These, however, are ac- 
tually unessential restrictions, and to illustrate this we 
have performed Brownian dynamics simulations for Dl 
of a non-truncated and rather long-ranged soft repul- 
sive potential (representative of highly charged colloids 
at low ionic strength), namely, the Yukawa potential 
u{r)/kBT = Kexp[-z{r/a - l)]/{r/a) with K = 554 
and z = 0.149. 

These data are compared in Fig. 1 with the corre- 
sponding data for the TLJ system (z/ = 6), much closer 
to the HS limit (1^ = 00, Ref. [HI), also shown in the 
figure. For the three systems we plot 1/D* as a func- 
tion of (j) (inset) and of the effective HS volume fraction 
(t^Hsin-i T, v), which is obtained not from with blip func- 
tion method [l3|, but from the iso-structurality condi- 
tion, written as S{kmax;n,T,iy) = SnsikmaxO-HS] 4>hs)- 
This condition requests that the hight of the main peak 
of the static structure factor of the "real" soft-sphere sys- 
tem, and of the effective hard-sphere system, coincide. As 
observed in the main figure, the data for D* of the three 
systems indeed collapse reasonably well when plotted as 
a function of (j)Hs{n,T^v). Let us also notice that the 
self-consistent generalized Langevin equation (SCGLE) 
theory of colloid dynamics (Ecs. (1), (2), and (5-8) of Ref 
17| . with kc = l.SSfcmaa:), Complemented with virtually 



exact liquid-theory approximations for S{k), provides an 
excellent first-principles quantitative description of these 
data without any adjustable parameter (solid lines of Fig. 

In the same figure we have also included some exper- 
imental data for the relative low-shear viscosity rj* = 
r]{n)/r]{n = 0) of several microgel soft-sphere solutions 
reported in Ref. [l2|- These data correspond to samples 
with different soft-sphere size and softness which, upon a 
linear rescaling of the concentration, collapse onto a mas- 
ter curve (Fig. 2(a) of Ref. [l^). As illustrated in Fig. 
1, such experimental master curve for rj*{(j)) coincides 
pretty well with our HS simulation data for l/D*{(j)Hs), 
thus indicating that these samples clearly belong to the 
hard-sphere dynamic universality class. One important 
question then refers to the conditions under which a sys- 
tem with arbitrary repulsive interaction will belong to 
this dynamic universality class. Our conjecture is that 
such condition is the existence of a distance of closest 
approach amim such that u^'^\r) ^ fc^T for r < amim 
so that particle-particle overlaps are highly unlikely or 
forbidden. Thus, according to this conjecture, systems 
with ultrasoft repulsive interactions with finite overlap 
potential energy u'^'^\Q) will not belong to this HS dy- 
namic universality class if fcsT w u'^'^\Qi). To illustrate 
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FIG. 1: Simulated data (solid symbols) and SCGLE theo- 
retical predictions (solid lines) of the normalized long-time 
self-diffusion coefficient D* {n,T,v), as a function of (inset) 
and as a function of (pHsi'n, T, v) (main figure) for the repul- 
sive Yukawa fiuid (squares) , the truncated 6-12 Lennard- Jones 
fluid (triangles) and the hard sphere fluid (circles, Ref. [ij]). 
The empty diamonds and the asterisks represent, respectively, 
the experimental master curve from Fig. 4 of Ref. 12] for 
the reduced viscosity (proportional to 1/D*) of microgel solu- 
tions, and the data for Dl{4i) of third- and fourth-generation 
dendrimer solutions in Fig. 3 of Ref. [l^ , scaled to collapse 
among themselves and with the HS master curve at low den- 
sities. 



this notion, in Fig. 1 we have also included the experi- 
mental data for D* of two low-generation dendrimer so- 
lutions, which cannot be adjusted by our HS universal 
curve, and whose structural pro per ties are best described 
by the Gaussian core potential 
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As illustrated in Fig. 1, the SCGLE theory accurately 
captures the scaling rules implied by Eqs. ^ and (H)). 
In reality, most of the work leading to the results pre- 
sented here was actually guided by the general qualitative 
predictions of this theory. For example, from the gen- 
eral mathematical structure of the theory, it is clear that 
the dynamic equivalence principle just illustrated implies 
a far more general density-temperature-softness scaling, 
best described in terms of the universal iso-dynamical 
surfaces in the (n, T, state space. These are defined as 
the loci of the points with the same dynamical proper- 
ties. According to Eqs. ^ and (jH), the iso-diffusivity 
surfaces labeled by the condition D*{n^T,v) = I?*, 
with T)* being a prescribed constant value, are also iso- 
T* and iso-(/>if5 surfaces. Referring to the TLJ model 
system, in Fig. [2] we plot the cuts along the plane 
^' = 6 of these iso-dynamical surfaces corresponding to 
V* = 10"\ 10"^ 10^3 an^i g.O (calculated, however, 
as \so-(f)HS lines such that 4>hs{4'tT* ,1^) = 0.494, 0.55, 
0.58, and 0.582, respectively). The correspondence be- 
tween the labels (j)Hs and V* was based, for the case 
of (j)HS — 0.494 (i-e^ V* « 10~^), on Lowen's dynamic 
freezing criterion |18| , and for the other iso-dynamical 
curves on the predictions of the SCGLE theory of col- 
loid dynamics, which locates the ideal glass transition at 
^HS = 0.582 [nl. 



FIG. 2: State space of the truncated 6-12 Lennard- Jones fluid, 
showing four iso-diffusivity lines, including the freezing tran- 
sition (dash-dotted line) and the ideal glass transition (solid 
line). In the inset we plot similar iso-diffusivity lines for the 
soft-sphere system with harmonic repulsive potential of Ref. 
[^, together with the simulation results for the Iso-Tq lines 
labeled = 10\ 10^, 10^, and 10* in Fig. 1 of the same ref- 



Let us mention that it has recently been discov- 
ered [l^l that the SCGLE theory of colloid dynam- 
ics can also describe the long-time dynamics of atomic 
systems, provided that the short-time self-diffusion co- 
efficient Z?" assumes its kinetic-theory value = 
{^/n/16ip)[a^/hyBT/M]. This suggests the manner in 
which the density-temperature-softness scaling and the 
iso-dynamical scenario just discussed, may be shared by 
atomic systems. Notice, however, that in this case be- 
comes state-dependent, and hence, in contrast with the 
Brownian case, we cannot expect that an iso-r* surface 
will also be an Iso-Tq, surface. To illustrate this, in the 
inset of Fig. [5] we compare a few SCGLE-predicted iso- 
T* lines for a soft-sphere fluid with harmonic repulsive 
potential u{r) = e(l — r/tr)^ for r < a, with the Iso-Tq 
lines determined by molecular dynamics simulations by 
Berthier and Witten [5] for this model system. In spite of 
the expected quantitative differences, the theoretical and 
simulated scenarios are, however, qualitatively identical. 

Fig. [1] provides the first illustration that the results for 
D'^g[(j)\ actually play the role of universal master curves 
for all moderately soft sphere systems, provided that the 
horizontal axis is labeled not by the volume fraction 
but by the effective volume fraction 0_f/5 (0, T* ,v). This 
scaling, however, can be expressed in alternative man- 
ners. Thus, we can plot D*fjg[(f>\ (or T^5[fccr;0]) not as 
a function of but as a function of a combination of 
such as [(t>Zcs{<f)\- The resulting curves will then also 
be master curves for all the soft sphere systems in the 
HS dynamic universality class when the horizontal axis 
is not the combination [(pZcsi't')] bnt the combination 

A very important meaning of the resulting master 
curves is suggested by the structure of the equation of 
state in Eq. ([T|), which can also be written in terms of 
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FIG. 3: Theoretical (solid line) and simulated (symbols) hard- 
sphere master cmve for D^si't'Hs) from Fig. [U but now 
plotted as a function of the combination [<j}HsZcs{(t>Hs)] = 
\l{T*)p*/T' (main panel) and T* / {p* \i{T*)) (inset). 



the dimensionless pressure p* = iKJ^py/Qt as 
pl{(l),T*)Xl{(j),T*) 

= (pHs[<P,T ,v)Zcs{(PHs{(P,T ,v)) 

(5) 

with Xu((j),T*) = fT//s((/), T*, z^)/(T. Thus, the combina- 
tion [4>HsZcs{4>Hs)] is also the combination [p*Xl/T*]. 
In Fig. |3] we replot the master curve Dj^glcfl of Fig. 
[1] but now as a function of this combination. Alterna- 
tively, we can use the combination [(l)Zcsi4')]~^ •> which 
changes the horizontal axis to [T*/p*X'l]. In the inset 
of Fig. |3] we replot the same results for D'^jglcj)] as a 
function of this combination. The resulting curves in 
Fig. El then constitute a prediction of the existence of 
master curves that describe the universal dependence of 
D*{(j),T* on pressure and temperature for this class 
of soft-sphere systems. The prediction of similar master 
curves for T*{ka;(l),T*,v) can be drawn from the deter- 
mination of the results for T^^gika; cj)). 

This predicted scaling has actually been corroborated 
by the master curve empirically discovered by Xu et al. 
[J| in their recent simulations performed on the TLJ soft 
sphere system with v = 6 and on other soft potentials. 
These authors found that the simulated results for the 
a-relaxation time Tq, of these soft-sphere systems, com- 
puted as a function of temperature T at fixed pressure 
p, or as a function of 1/p at fixed T (Figs, l.a and l.b 
of Ref. [il , respectively) all collapse onto a master curve 
whenplotted as a function of the ratio T/p (Figs. 2 and 
3 of y). The simplest manner to relate their empirical 
scaling with our predictions in Fig. |3] above, is to scale 
the raw data for Tq in their Fig. 1 as r* = k'^D^Ta, 
with £)° = {s/TT/l&4>)[(j^Jk^TjM], and to plot them as a 
function of [T* /p*Xl{T*)]. As a result, the various data 
of Fig. 1 of Xu et al. ^4] collapse onto the master curve 
shown here in Fig. |4j predicted by the dynamic equiv- 
alence discussed in the present work, and in agreement 
with the scaling discovered by Xu et al. This required 
to express (j) in this expression for in terms of T* 
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FIG. 4: Molecular dynamics simulations reported in Figs, la 
and lb of Xu et al. 01 for Tc of soft-sphere liquids, plot- 
ted here as the master curve for T*{T*,p*) as a function 
of the combination [<j)Hs{<l},T'' ,v)Zcs{<I}hs{<I},T* ,y))]^^ = 
\T* / \^{T*)p*\. In the inset we show the simulation data, ob- 
tained by Xu et.ai, for the p* /T* ratio (squares), compared 
with Eq. (O (solid line). 



and p* , but this is easily achieved using the equation of 
state in Eq. ([5]); the inset of Fig. |4] compares our theo- 
retical equation of state (Eq. jp) with the correspond- 
ing simulation data of Ref. Let us finally notice 
that for the TLJ model, in the low-temperature regime 
X^((t>. T*) may be approximated by the blip-function re- 
suh [i, [l^, which yields Xl(T*) « 1 - (30F/2:y)VT^. 
Thus, at the temperatures employed in the simulations 
of Xu et al. (T* < lO'^), the combination [T* /p* Xl(T*)] 
is essentially the temperature-to-pressure ratio [T*/p*], 
employed in Ref. 

In summary, we have illustrated the accuracy of the 
density-temperature-softness scaling of the dynamics of 
soft-sphere liquids, and shown that it provides a simple 
and useful conceptual tool to understand, within a unified 
framework, the phenomenology of both, thermally-driven 
molecular glass formers and density-driven hard-sphere- 
like colloidal liquids. There are, of course, important 
pending assignments, such as the characterization of the 
non-equilibrium dynamics at, and beyond, the glass tran- 
sition in this (HS) universality class, or the identification 
of other dynamic universality classes, particularly those 
influenced by the presence of attractive interactions. We 
expect that the results discussed in this work, together 
with the non-equilibrium extension of the SCGLE theory 
[20j . will facilitate the progress in these directions. 
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